Physics of the Pseudogap in 8-site Cluster Dynamical Mean Field Theory: 
photoemission, Raman scattering, in-plane and c-axis conductivity 
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Cluster dynamical mean field and maximum entropy analytical continuation methods are used to 
obtain theoretical estimates for the many-body density of states, electron self-energy, in-plane and 
c-axis optical conductivity and the Big and B2g Raman scattering spectra of the two dimensional 
square lattice Hubbard model at intermediate interaction strengths and carrier concentrations near 
half filling. The calculations are based on an 8-site cluster approximation which gives access to both 
the zone-diagonal and zone-face portions of the Fermi surface. At low dopings the zone-face regions 
exhibit a 'pseudogap', a suppression of the many-body density of states for energies near the Fermi 
surface. The pseudogap magnitude is largest near half filling and decreases smoothly with doping, 
but as temperature is increased the gap fills in rather than closes. The calculated response functions 
bear a strong qualitative resemblance to data taken in the pseudogap regime of high- Tc cuprates, 
strongly suggesting that the intermediate coupling Hubbard model accounts for much of the exotic 
behavior observed in high-Tc materials. 

PACS numbers: 71.10.Fd, 74.72.-h, 71.27.-|-a,71.30.-|-h 



I 

o 
o 



> 

^' 

O 
O 



I. INTRODUCTION 

The 'pseudogap', a temperature and carrier concen- 
tration dependent suppression of the many-body den- 
sity of states of hole-doped high temperature copper- 
oxide superconductors which is visible at temperatures 
well above the superconducting transition temperature, 
is one of the enduring mysteries of the field. The pseu- 
dogap was first inferred from measurements of spin- 
lattice relaxation timesi and Knight shiftsj^ additional 
evidence rapidly accumulated from measurements of in- 
plane resistivity)^ photoemission spectra?^ the inter- 
plane (c-axis) conductivityj^ii Raman scattering^ and 
tunnelingi^ and measurements of the electronic density 
of states. On the other hand, the pseudogap docs not 
lead to a significant suppression of low frequency spectral 
weight in the in-plane conductivityjii although structure 
in the conductivity and in scattering rates inferred from 
the conductivity has been attributed to the pseudogapii^ 
The gap is most pronounced at low temperatures and low 
dopings. The data suggest that the gap magnitude de- 
creases as doping increases, whereas with increasing tem- 
perature the gap magnitude docs not decrease. Rather, 
the gap 'fills in' as more and more states appear in the 
mid-gap region. 

Photoemission measurements^i^ii^ indicate that the 
pseudogap is largest near the zone corner (0, 7r)/(7r, 0) 
regions and vanishes at the zone-diagonal (±7r/2, ±7r/2) 
regions. This momentum-space structure of the pseu- 
dogap is one aspect of the more general phenomenon 
of momentum-space differentiation that characterizes the 
doping-dependent metal insulator transition in the high- 
Tc materials. 

Schmalian, Pines and Stokjovio^^ argued that one 
should distinguish between a 'weak' and 'strong' pseudo- 
gap. The weak pseudogap terminology refers to relatively 
low energy phenomena affecting states within a few tens 



of meV of the Fermi surface while the 'strong pseudogap' 
is a suppression of density of states which may exist over 
a relatively wide energy range. 

The physical interpretation of the pseudogap remains 
controversial. One possibility is that it is associated with 
an actual thermodynamic transition to a phase with a 
definite long ranged order. Possibilities which have been 
proposed include magnetic order, perhaps of 'stripe' or 
'nematic' density wave formic and an orbital current 
phaseii^ Experimental evidence has appeared support- 
ing each of these possibilitie a^^d^ but the interpretations 
remain controversialji^ 

An alternative possibility is that the pseudogap is a 
consequence of long but not infinite ranged order of spin 
density wave ;^^'^"" — superconducting^^— or RVI3^"— 
type. One loop>^ and more sophisticate d^^i^^'^^ cal- 
culations provide a relation between quasi-long-ranged 
order and pseudogaps, while slave-boson-based mean 
field methods^^"— provide a different mechanism for 
(pseudo)gap formation. However, these methods are 
based on uncontrolled analytical approximations. Mod- 
els of quasi one dimensional 'ladder' compounds can be 
studied in a controlled manner and are known to ex- 
hibit gapped phases with no long ranged order )^ however 
despite intriguing qualitative similarities the relation of 
'ladder' calculations to the two dimensional physics rel- 
evant to the cuprates remains unclear. There is a clear 
need for studies based on methods which are applicable 
in two dimensions and at intermediate to strong cou- 
plings and which are not based on a particular assump- 
tion about the type of relevant correlations. 

One such approach is provided by 'cluster' dynami- 
cal mean field methodsi^ These techniques are based on 
a coarse discretization of the momentum dependence of 
the electron self-energy but permit a numerically unbi- 
ased solution of the resulting model. Important early 
work showed that the cluster dynamical mean field ap- 
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proximation produced features reminiscent of the pseu- 
dogap including suppression of the density of states in 
the (0, tt) region of the zone^ and 'Fermi arcs' which 
are at least qualitatively consistent with photoemission 
experiments)^^— Subsequently many cluster dynamical 
mean field theory (DMFT) studies of the pseudogap have 
appeared.—"— 

Studies to date are mainly of two sorts: large-cluster 
studies of one doping and temperature and more com- 
prehensive studies of smaller (two and four site) clusters. 
An important large cluster study is the work of Macridin 
et al.^ which considered a 16 site cluster at t/ = 8t at 
a doping of 5% and an inverse temperature /3 = 8/i and 
showed from analytical continuation of the electron spec- 
tral function that a (0,7r) pseudogap existed at this dop- 
ing. The smaller-cluster studies have been based on two 
and four site clusters for which the computational bur- 
den is much less, permitting systematic examinations of 
a wider range of parameter space . •^'^I'^^i'^^i'^^^^^i"^^' '— How- 
ever these smaller clusters do not allow direct compari- 
son of zone-diagonal and zone-face regions of momentum 
space. In addition, many of the more systematic small 
cluster based analyses employed interpolation schemes 
to construct a representation of the electron self-energy 
throughout the Brillouin zone. The reconstructions are 
based on data corresponding to the (0, 0), (tt, tt) and, for 
the four-site clusters, (0, 7r)/(7r, 0) momentum points^""— 
and the physical interpretation is problematic especially 
because information related to the physically important 
(±7r/2, ±7r/2) momentum space regions is inferred from 
measurements at (0, 0), (tt, tt) and (0, tt) / (tt, 0) where the 
physics is presumably different. However, an interesting 
recent paper by Ferrero et al.—i^ introduced a two-site 
cluster with a new momentum space partitioning which 
clearly separated nodal and antinodal regions. 

The work reported here builds on recent studies of the 
8-site cluster— 1^ shown in Fig. [T] The larger size of this 
cluster provides a more refined momentum resolution and 
in particular gives independent access to the zone-face 
and zone-diagonal regions of the Fermi surface. However 
the size is small enough that studies of wide regions of 
parameter space are computationally feasible. Previous 
work has revealed that the doping-driven Mott transition 
is momentum-space-selective, with a gap opening first in 
the zone corner (0, tt)/{tt, 0) regions of the Brillouin zone 
while the zone-diagonal (±7r/2, ±7r/2) regions remain un- 
gapped until the carrier concentration reaches half filling. 
The previous work focused on quantities defined directly 
on the Matsubara axis. In this paper we use analytical 
continuation techniques to examine the consequences of 
the momentum-space-selective transitions for observables 
including the electron spectral function and self-energy, 
the interplane and in-plane conductivity, and the Raman 
scattering intensity. Our results strongly suggest that 
even in the absence of long ranged or quasi-long-ranged 
order, the Hubbard model at intermediate couplings con- 
tains the essential physics of the 'strong pseudogap'. 

The rest of this paper is organized as follows. In Sec- 



tion |TT] we present the model and the dynamical mean 
field and analytical continuation methods we use to solve 
it. In Section ITTll we show results for the electron spectral 
function and in Section IIVI the interplane conductivity. 
Section |V] presents our results for the Raman scattering 
intensity and Section IVII the in-plane conductivity. Sec- 
tion IVIII shows the electron self-energy in the different 
momentum sectors, confirming the conjecture that the 
gap arises from an orbitally selective Mott transition and 
demonstrating that the model reproduces key aspects of 
the momentum selectivity in the approach to the Mott 
transition. Section [Villi is a summary and conclusion. 

II. MODEL AND METHODS 

We study the Hubbard model on a two dimensional 
lattice. The model is conveniently written in a mixed 
momentum-space/position-space representation as 

H = '^Skcl^Cka + U ^n^t^-ii - ^i^^lMa. (1) 

ka i i(7 

We take, as a reasonable representation of the band 
structure of high temperature superconductors, 

Ek = — 2t (cos kx + cos ky) — At' cos kx cos ky (2) 

with t' = — 0.15i. For comparison to data we note that 
a value generally accepted for high- Tc superconductors 
is t ss 0.35 eV^ while values of t'/t from -0.05 to -0.3 
have been reported for different materials. 

To solve the model we use the 'dynamical cluster ap- 
proximation' (DC A) 1^1^ The method is based on tiling 
the Brillouin zone into N equal-area non-overlapping tiles 
and approximating the electron self-energy T,{k,ijj) as a 
piecewise constant function which may take different val- 
ues in the different tiles. Labeling the tiles by center 
momentum K we have 

S](fc,u;) = SA'(a;); k € K. (3) 

The results we present were obtained using the 8- 
site momentum space partitioning shown in Fig. [T] and 
the 'continuous-time auxiliary field' (CT-AUX) numeri- 
cal method^ as discussed in more detail in Ref. and 
[sol . Because the model is solved on the imaginary axis, 
an analytical continuation procedure is required to ob- 
tain real frequency information. Following Ref. [s^ we 
continue the electron self-energies using the Maximum 
Entropy technique^ and the 'L-curve' method.— The 
covariancc matrix of the self-energies is approximately 
diagonal and the continuation of the obtained real fre- 
quency spectra back to the imaginary axis is in good 
agreement with the original data. Although uncertain- 
ties exist in the analytical continuation, our experience 
is that the near Fermi-surface structures are reliable. 

The 8-site computation suffers from a fermionic sign 
problem which becomes worse with decreasing tempera- 
ture and increasing doping and interaction strength. The 
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3. Quasiparticle contribution to Raman scattering in- 



FIG. 1: Brillouin zone partitioning associated witli tiie 8-site 
cluster. 



need to scan a range of dopings and to obtain data of 
the precision required for reliable analytical continua- 
tions limited us to temperatures > i/20 (corresponding 
roughly to 200 K for t 0.35 eV) and U = It. The 
[/-value is smaller than the values U ^ 9t believed^ to 
describe high temperature superconductors, but is large 
enough that (within the approximation we use) the half 
filled state is a Mott insulator and hole-doping leads to 
a momentum-space-selective Mott transition. 

Dynamical mean field methods involve a drastic sim- 
plification of the momentum dependence of the electron 
self-energy. As Eq. ^ shows, the DCA method produces 
a piecewise constant self-energy which may be viewed 
as a discretization of the continuous momentum depen- 
dence of the exact solution. In many previous dynami- 
cal mean field studies of the pseudogap an interpolation 
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was used to construct a self-energy with a 
continuous momentum dependence, which was then used 
to produce figures to be compared to experiment. We 
prefer to avoid interpolations and instead work with ana- 
lytical continuations of directly calculated quantities. We 
present results for 

1. Sector-averaged electron spectral function 



Nr. 
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(4) 



with Nc = 8 in our 8-site cluster. 
2. Interplane or c-axis conductivity 



f d k f dcj 



t^{k)lmGK{k,uj + n)t^{k)lmGK{k,uj) (5) 

with t±{k) = to(cosfc^ — cosfc.y)^, GK{k,uj) = 

[uj ~ Ek + fJ- — Sa'('^)] and /(w) the Fermi-Dirac 
distribution. We used to = 0.15 eV w 0.434^ in 
the calculation. 



tensity 
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Two geometries are of interest: Big, where 
matrix element is jab = ^{d^Sk/dk^ — d'^Sk/dky) 
and i?2g where jab = d'^Sk/ dkxdky. The Big 
geometry highlights the zone-face (0, tt) / (tt, 0) re- 
gions while the B2g geometry highlights the zone- 
diagonals (±7r/2, ±7r/2). This matrix element is 
appropriate for non-resonant Raman scattering and 
is the simplest one which is consistent with the sym- 
metry. In practice incident laser frequencies are of- 
ten chosen to take advantage of resonant enhance- 
ments arising from other degrees of freedom in the 
solid. These will change the absolute and relative 
magnitudes but not the symmetries of the vertices. 

4. Quasiparticle contribution to in-plane optical con- 
ductivity 



d^k f duj f{io) - f{LO + n) 



VxlTnGKik, uj -t- ^T}VxhnGK{k, ui) 



(7) 



where 



dek/dk-x- An approximate vertex 



correctioiias was also incorporated. 

5. Real and imaginary parts of sector-dependent elec- 
tron self-energy EA-(aj). 

Here denotes an integral over momenta lying in sec- 
tor K . Note that for the Raman scattering and in-plane 
optical conductivities (unlike for the other quantities we 
have considered) a vertex correction contribution (which 
we have only partially calculated) is present. The full 
vertex correction calculation is currently in progress us- 
ing methods outlined in our previous work^ but based 
on these results we expect the low frequency conductiv- 
ity of primary interest here to have only a small vertex 
correction. 



III. ELECTRON SPECTRAL FUNCTION 

Fig. [2] shows the electron spectral function for the sec- 
tor containing the (0, tt) momentum (labeled G in Fig.[T]) 
calculated at 0.05 hole doping for several temperatures. 
A 'pseudogap' (reduction in density of states) is visible 
in the low energy region. 

We define the pseudogap magnitude Epc = ^Apo as 
the peak to peak separation (for the x = 0.05 case shown 
in Fig.^EpG w 0.8t w 0.3 eV). The reduction in den- 
sity of states is largest at the lowest temperature and for 
frequencies near w = 0. It appears that at this doping 
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FIG. 2: Many-body density of states A{Ltj) = — ^lmG'(tj) av- 
eraged over sector C containing the (0, vr) point, calculated 
using 8-site DC A approximation at hole doping x = 0.05 with 
U = 7t and inverse temperatures /3 = l/T indicated. 



the low frequency density of states vanishes as T — !■ 0. 
As temperature is raised the gap 'fills in': the density 
of states inside the gap increases but the gap magnitude 
does not change appreciably. For temperatures greater 
than about 0.2t the gap is no longer visible at this doping. 

The upper panel of Fig. [3] shows the (0, 7r)-sector spec- 
tral function calculated at several different dopings. A 
decrease in gap magnitude with increasing doping is evi- 
dent. For dopings larger than x = 0.11 a gap is not visible 
at the temperatures T > t/20 accessible to us although 
a weak feature in the x = 0.13 curve suggests that the 
gap is still present. However, certainly at x = 0.11 and 
perhaps at a; = 0.13 the gap magnitude (as defined by 
the peak-to-peak distance in the spectral function) is not 
small. We therefore suspect that at least a reduction of 
density of states would be observed at higher dopings if 
we were able to perform the calculations at lower tem- 
peratures. 

The lower panel of Fig. [3] shows the (7r/2, 7r/2)-sector 
spectral function at the same dopings. At the smallest 
doping a weak suppression of low frequency density of 
states is evident but for most dopings this sector remains 
ungapped. 



IV. INTERPLANE CONDUCTIVITY 

An important early indication of the presence of a 
charge-pseudogap was provided by measurements of the 
frequency dependence of the interplane conductivity)^ As 
can be seen from Eq. ([S]), in high- Tc materials the matrix 
elements relevant to the interplane conductivity highlight 
the zone-face regions where the electron spectral function 
exhibits a gap (see upper panel in Fig. 

Fig. m shows the calculated temperature and doping 
dependence of the interplane conductivity. The pseudo- 
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FIG. 3: Many-body density of states A{lj) = -ilmG(aj) 
averaged over sector C containing the (0, tt) point (upper 
panel) and sector B containing the (7r/2,7r/2) point (lower 
panel), calculated using the 8-site DCA approximation at 
hole dopings indicated with U — 7t and inverse temperature 
/3 = 20/t 200 K. 



gap is visible as a temperature and doping dependent 
suppression of the low frequency interplane conductivity. 
The interplane conductivity is suppressed over a rela- 
tively wide frequency range; the suppression increases as 
the doping or temperature decreases, and the gap fills 
in but does not close as temperature is increased. The 
calculations also reveal a weak maximum in the conduc- 
tivity at an energy just above the pseudogap. A some- 
what broader version of this feature was observed by Yu 
et al^ It is possible that the relative sharpness of the 
feature is an artifact related to our coarse-graining of 
momentum space, which might arise because the DCA 
approximation necessarily produces a gap that is piece- 
wise continuous; and as is known from the familiar case 
of s-wave BCS superconductivity a momentum indepen- 
dent gap produces a peak. The results are reasonably 
consistent with experiment i^i^ i^^'^^ Ref. [s^ reports a high 
energy pseudogap of a magnitude consistent with what 
is found here. It is important to note that in the widely 
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FIG. 4: Upper panel: temperature dependence of interplane 
conductivity (Tc(n) calculated from Eq. ([SJ for hole doping 
X = 0.05 at inverse temperatures /3 = l/T indicated. Lower 
panel (note change of vertical scale): doping dependence 
of interplane conductivity calculated at inverse temperature 
P = 20/f ^ 200 K. At hole doping x = 0.16, the low fre- 
quency interplane conductivity is found to be approximately 
of Lorentzian form with dc value 0.65. 

studied YBa2Cu306-(-a; material the interplay of strong 
local field effects (arising from the bilayer structure) and 
phonon effects produce complicated structures in the low 
frequency conductivity which arc not represented in the 
present calculationi^"— 

Conductivities may be characterized by 'spectral 
weight', the integrated area in some frequency range. 
The total spectral weight obeys an 'f-sum' rule, which 
for the model studied here is 




f{uj)t\{k)lm[G\k,uj)]. 

We have verified that the spectral weight obtained from 
integration of the conductivity equals the spectral weight 
obtained from an evaluation of Eq. ([S]) . The curves shown 
in Fig. |4] imply a temperature and doping-dependent de- 
crease in low frequency spectral weight. We find that in 
the lower doping regions where the pseudogap is present, 
the spectral weight which is lost at low frequencies due 
to the formation of the gap is not fully restored at 
higher frequencies; thus the conductivity spectral weight 
decreases as the temperature or doping is decreased, 
whereas at higher doping the c-axis conductivity spec- 
tral weight increases as temperature is decreased. The 
calculated interplane conductivity is spread over a wide 
frequency range so the pseudogap-induced decrease of 
spectral weight is small relative to the total weight. The 
doping and temperature dependence of the c-axis spec- 
tral weight is given in Table ID 

Ref. |48| reports interplane conductivity results obtained 
using a two-site cluster, a particular choice of self-energy 



periodization and a Fade continuation to compute Uc for 
U/t = 10 and t' /t = —0.3. The results are very similar 
to those shown here but with slightly larger gaps and a 
much greater suppression of CTc at sub gap frequencies. 
Note that a factor of Trig/ (8i^) is required to convert the 
results of Ref. to our units. 



X pt ^ 10 I3t = 15 pt = 20 
0.05 0.111 0)5 0.102 
0.07 0.125 0.119 0.116 
0.08 0.142 0.138 0.134 
0.16 0.207 0.223 0.231 



TABLE I: Integrated spectral weight of interplane optical con- 
ductivity obtained from Eq. at different temperatures for 
several hole dopings. 

To conclude this section we note that the pseudogap 
effects are enhanced by the structure of the matrix el- 
ement (which is the one widely used in the theoretical 
literature) which enhances the contribution of sector C 
relative to sector B. If a fc-indepcndent matrix element 
was used, the pseudogap effects would be much less pro- 
nounced. 



V. RAMAN SCATTERING INTENSITY 

Raman scattering has provided important information 
about the cuprates.— It is a photon- in/photon-out pro- 
cess, and the polarizations of the electric field vectors 
of the incident and emitted photons may be adjusted to 
highlight transitions in different regions of momentum 
space. We consider here two scattering geometries: Big 
which highlights the anti-nodal region, and B2g where the 
matrix element highlights the nodal region. The compu- 
tation of the Raman intensity involves a vertex correc- 
tion, which we have not computed (note that unlike the 
case of the in-plane conductivity^ the vertex correction 
does not involve a contribution from the momentum- 
space discontinuity of the self-energy). Thus what is 
shown is the quasiparticle contribution only. For this 
reason, and also because we do not have control over the 
high frequency part of the spectrum, we do not undertake 
the sum rule analysis introduced by de Medici et al^ in 
their single site DMFT study of the Raman intensity in 
the strong correlation limit. 

Fig. [5] shows the temperature dependence of the cal- 
culated Raman spectra in the two scattering geometries 
at the low doping a; = 0.05 (upper panels) and high dop- 
ing X — 0.16 (lower panels). We see that at the higher 
dopings the functional forms and temperature dependen- 
cies in the two scattering channels are similar, while at 
the lower doping they are rather different. At the higher 
doping the Raman response in each channel rises linearly 
to a weak maximum and then approximately saturates; 
the initial slope increases as T decreases in both chan- 
nels. It is also interesting to note that the calculation 
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FIG. 5: Temperature dependence of Raman scattering inten- 
sity x"(^) calculated from Eq. for Big and B-2g scattering 
channels in 8-site DCA at hole doping 0.05 (upper panels) 
and 0.16 (lower panels) and U — It. 



reproduces the roughly frequency-independent high fre- 
quency behavior, which had previously been argued to 
be evidence for novel 'marginal Fermi liquid' physics not 
contained in the Hubbard model<^ However, at the lower 
dopings the low frequency Big response is progressively 
suppressed as T decreases, whereas the B2g response has 
negligible T dependence. 

The key features of the results presented in Fig. [SJ 
namely an increase with increasing doping in the tem- 
perature dependence of the low frequency i?2g scattering 
intensity and a change in sign of the T dependence for the 
Big intensity along with the presence of a weak maximum 
at an energy which decreases with increasing doping arc 
qualitatively consistent with data (see e.g. Fig. 42 and 
43 of Ref.lsi). 

However, there are important quantitative differences. 
The calculated ratio of Big to B2g is approximately 15, 
whereas the experimental ratio is much closer to one. 
This difference presumably arises because the experi- 
ments employ a resonantly enhanced matrix element. 



More importantly, the calculation places the maximum 
(or scale at which the response saturates) at a higher en- 
ergy than the data in Refs. [6^ and [g^- For example at 
X = 0.16 the calculated saturation point for the Big spec- 
tra (using t = 0.35 eV and placing the saturation maxima 
at 0.33t) at /3t = 20 is w w 925 cm^^ while in the Bag 
spectra our results saturate at 700 cm~^ (placing the 
saturation point at 0.25i) whereas the experimental data 
have not quite saturated. Further, observed B2g Raman 
spectra at low dopings have more T dependence than is 
found in our calculations. The onset coincides in tem- 
perature with other measures of the pseudogap, so we 
associate the phenomenon to the pseudogap. 

Figs. [5] and [7] show the doping dependence of the Ra- 
man intensity at inverse temperature pt = 20 corre- 
sponding to T w 200 K. A strong increase in initial 
slope is evident in both channels. The rapid steepening 
of the initial slope in the i?2g channel is a consequence 
of the emergence of coherent quasiparticles in the zone- 
diagonal sector (see Sec. IVIip . The change in the Big 
channel arises from the doping dependence of the pseu- 
dogap. The calculated behavior of the Raman intensity, 
including the difference in doping and temperature de- 
pendence between the two sectors and in particular the 
doping dependent suppression of the Big Raman spectra 
over a wide frequency range, is in reasonable agreement 
with measurements . However, our calculated i?2g 
spectra exhibit a strong doping dependence which is not 
observed in recent experimentsi^ 
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FIG. 6: Raman scattering intensity x"(^) in the -Big geom- 
etry at different hole dopings as indicated with U = It and 
inverse temperature j3 = 2Q/t ~ 200 K. 



VI. IN-PLANE CONDUCTIVITY 

The optical conductivity uljl) of the high- Tc cuprates 
has been an enduring mystery. The salient features 
of the data are a 'Drude' peak centered at zero fre- 
quency, with a strongly temperature dependent spectral 
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FIG. 7: Raman scattering intensity in the B2g geom- 

etry at different hole dopings as indicated with U = 7t and 
inverse temperature /? — 20/t ~ 200 K. 



FIG. 9: Contribution of different momentum sectors to cal- 
culated in-plane optical conductivity for U — 7t at doping 
X = 0.05 and inverse temperature /3 = 20/t ~ 200 K. 
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FIG. 8: Main panel: in-plane optical conductivity a{Q,) cal- 
culated from Eq. ^ at U = 7t and inverse temperature 
/3 = 20/t ~ 200 K for hole dopings indicated. Inset: tem- 
perature dependence of the low frequency in-plane optical 
conductivity for hole doping x — 0.16. 



tivity. Our results exhibit strong similarities to exper- 
imental data, including, at low frequencies, a growth 
of the 'Drude' peak with doping, and a weakly fre- 
quency and doping-dependent 'mid-infrared' conductiv- 
ity with a magnitude comparable to the measured value 
~ 5 — 700 n^^cm~^ (note that to convert to the unit 
of fl~^cm~^ commonly used in experiments, our results 
have to be multiplied by a factor of about 5 x 10"^). The 
inset shows that at higher dopings the 'Drude' peak grows 
noticeably and sharpens slightly as temperature is de- 
creased. 



X fit ^10 fit = 15 pt = 20 
0.05 0.252 0244 0.240 
0.08 0.3141 0.3145 0.313 
0.11 0.328 0.343 0.345 
0.16 0.383 0.390 0.392 



TABLE II: Integrated spectral weight of in-plane optical con- 
ductivity obtained by integrating calculated curves multiplied 
by a factor of 2/7r up to $1 = 3t « 1 eV at different tempera- 
tures for several hole dopings. 



weight and half width, and a broad higher frequency 
continuum ] ^ i , 12^69 j^^^^ strong doping dependence of the 
'Drude' peak has been taken as evidence of strong 'Mott' 
correlations, while the broad higher frequency continuum 
has been interpreted in terms of scattering from spin 
fluctuations,^i^-I^ 

The computation of the in-plane conductivity involves 
a vertex correction^ which has not been fully calculated. 
Here we only include the vertex correction arising from 
the self-energy discontinuities in momentum-space, while 
the contribution of additional dependence of the self- 
energy on the vector potential A which arises from the 
change in mean-field function is neglected. 

The main panel of Fig. [S] presents the calculated dop- 
ing dependence of the real part of the in-plane conduc- 



A particularly interesting issue is the lack of a clearcut 
effect of the pseudogap on the conductivity. A natu- 
ral conjecture is that the conductivity is dominated by 
states near the zone-diagonal, which are not sensitive to 
the pseudogap. Fig. |9] presents a decomposition of the 
contributions of the different sectors to the measured con- 
ductivity which supports this conjecture. A pseudogap 
(the feature at w t) is present in the contribution of 
sector C, but the contribution of this sector to the mea- 
sured conductivity is relatively small, so the feature is not 
evident in the full conductivity. The pseudogap feature is 
more pronounced in 4-site cluster calculation s^^' '''^'^^ be- 
cause the geometry of the 4-site cluster is such that the 
low frequency behavior is dominated by the (0, 7r)/(7r, 0) 
sectors so that it docs not capture the physics of the 
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nodal quasi-particles. In the present 8-site calculation 
the nodal quasiparticle physics is represented by sector 
B, which is seen to give the dominant contribution to the 
conductivity. 

The dopingiii^i^i^ and temperatur o^^'^^ dependence 
of the optical spectral weight K = J ^^a{uj) has been 
the subject of discussion in the literature. One ambigu- 
ity is the range over which the conductivity is to be inte- 
grated: setting the range too high leads to the inclusion 
of interband transitions which are believed to be irrele- 
vant to the physics of high- Tc materials while setting the 
range too low may mean that changes in the width of a 
low frequency peak may be mistaken for changes in area. 
Recent papers suggest that an upper cutoff of ^ 1 eV is a 
reasonable compromise value A marked increase of 
spectral weight occurs as doping is increased (see Ref. [13 
for a summary of the data) . The experimental consensus 
is that at all dopings the spectral weight increases weakly 
as temperature is decreased, and that the temperature 
dependence changes markedly as the temperature is de- 
creased below the superconducting transition tempera- 
ture. The doping dependence of the calculated spectral 
weight over the range f2 < 3t ^ 1 eV (shown in Table [ll]) 
is reasonably consistent with the data although we find 
that at lower dopings x = 0.08, 0.05 the temperature 
dependence flattens. A weak temperature dependence is 
seen in the spectral weights integrated up to 3t « 1 eV. 
The sign of the temperature dependence changes with 
doping: increasing as T is decreased at high doping and 
decreasing as T decreases at low doping. It is tempt- 
ing to relate this change in temperature dependence to a 
change in physics from Fermi-liquid-like at high doping 
to pseudogapped at low doping but the small magnitude 
of the effect and the possibility of temperature dependent 
changes in the functional form of the conductivity make 
an interpretation unclear. Calculations of a 4-site cluster 
approximation to the t — J model^i (which is believed 
to reflect the low energy physics of the Hubbard model 
for U > 12t, larger than what we studied here) found a 
similar doping and temperature dependence. 

The conductivity is sometimes expressed in terms of a 
frequency-dependent optical scattering rate Fopt related 
to the complex conductivity a = ai (w) -|- ia2 (w) via 



1.5 



opt 



KRe- 



(9) 



where cti is calculated from Eq. ([7]) and a2 is obtained 
by the Kramcrs-Kronig transformation of di . In under- 
doped cuprates at high temperature, Fopt is large and 
temperature dependent over a wide frequency range. As 
the temperature is decreased to the pseudogap scale the 
high frequency part loses its temperature dependence 
while at lower frequencies a temperature dependent sup- 
pression of Fopt appearsi^S Fig. [TU] shows that this be- 
havior is also found in our calculations. 
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FIG. 10: Optical scattering rate Fopt computed from Eq. Q 
for hole doping x = 0.11 and inverse temperatures /3 = 1/T 
indicated. 



VII. ELECTRON SELF-ENERGY 

Fig.[TT]and[T^comparc the temperature dependence of 
the imaginary part of the electron self-energy computed 
at doping x = 0.08 and x = 0.16 for the two momentum 
sectors containing the Fermi surface. At doping x = 0.08 
(see Fig. [TT|). we see that the (0, tt) sector self-energy is 
characterized by a pole located near cj = whereas no 
pole appears in the self-energy of the sector containing 
the zone-diagonal part. A near zero-energy pole in the 
self-energy is a characteristic of a Mott insulating state, 
confirming that the gap opening transition is indeed a 
sector selective Mott transition. As the temperature is 
decreased the pole grows in strength. The pole position 
at w w leads to the approximate particle-hole symme- 
try of the spectra. We also observe that at this doping 
the sector B (zone-diagonal) self-energy has only a weak 
temperature dependence at the temperatures accessible 
to us. 

At doping X = 0.16 (see Fig. [T^]), the self-energies in 
both sectors decrease with temperature and have a min- 
imum centered at w w 0. In sector B, the self-energy has 
Fermi-liquid-like behavior which decreases rapidly and 
roughly linear in T. In sector C, the self-energy decreases 
more slowly. The difference in magnitude and doping 
dependence indicates that this doping regime is charac- 
terized by a large variation in scattering rate around the 
Fermi surface. 

Fig. [13] shows the doping dependence of the imagi- 
nary part of the sector C self-energy. A pole, gaining in 
strength as doping is decreased and centered at approx- 
imately zero frequency, is clearly visible for x < 0.11, 
while for x = 0.13 and 0.16 there is no pole, only a 
weak modulation indicating a non-Fermi-liquid scatter- 
ing rate. Whether this modulation would evolve into a 
pole as T ^ is an interesting open question. Liebsch et 
al^ analysed the self-energy pole structure in the (0,7r) 
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FIG. 11: Imaginary part of electron self-energy computed at 
different temperatures for hole doping x — 0.08 in (0, 7r)-sector 
(upper panel) and (7r/2, 7r/2)-sector (lower panel). 



sector of a 4-site CDMFT study, finding a similar doping 
dependence of the pole strength. They reported a strong 
dependence of the pole position on doping; this variation 
is not found in the 8-site cluster studied here. 

An alternative characterization of electronic behavior 



1 - 



is the quasiparticle residue Z = 

Because within each sector the self-energy is momentum- 
independent, Z gives the renormalization of the Fermi 
velocity as v* = Zv. This renormalization has physical 
significance if the self-energy is Fermi-liquid-like, mean- 
ing that the imaginary part is not too large and the real 
part is linear in frequency over a reasonable range about 
oj = 0. We determine the boundaries of the Fcrmi- 
liquid regime by first observing that the real part of 
the self-energy ReE is linear in frequency over the range 
—cul < w < ujh, and then comparing the magnitude of 
the imaginary part of the self-energy at zero frequency 
to the change of w — (ReE(a;) — ReE(O)) over the linear 
range. If the change {ujh — RcT,{ujh)) — {ujl — ReS(a;L)) 



3, 



E 




FIG. 12: Imaginary part of electron self-energy computed at 
different temperatures for hole doping around 0.16 in (0, vr)- 
sector (upper panel) and (7r/2, 7r/2)-sector (lower panel). 
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FIG. 13: Imaginary part of electron self-energy computed for 
sector C containing the (0, tt) point at inverse temperature 
j3 — 20/t 200 K for dopings indicated. (For the real part 
see Fig.fTSl) 
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is larger than — 2ImS(aj = 0) we identify the regime as 
Fermi-hquid-hke. For an ihustration of the determination 
of Fermi- hquid behavior see appendix El 

This condition is reasonably well satisfied for sector 
B for dopings x > 0.08 (and marginally satisfied for 
X = 0.08). Similarly sector C is fomid to be Fermi-liquid- 
like for dopings x = 0.18 and greater, but for x = 0.06 
the self-energy in both sectors is far from Fcrmi-liquid- 
like and the quantity Z cannot be interpreted as a "quasi- 
particle weight" . 

The solid points in the upper panel of Fig. [T4l show the 
value of Z for the sector B containing the zone-diagonal 
point (7r/2,7r/2) and the sector C containing the zone- 
face point (0, tt) for dopings for which the sectors are 
Fermi-liquid-like. The open symbols show the mathe- 
matically defined values of Z in the regime where it has 
no physical meaning because the regime is not Fcrmi- 
liquid-likc. For dopings in the Fermi-liquid regime the 
Z in sector B is linear in x but extrapolates to a small 
non-zero value at a; = 0. This is approximately but not 
exactly the behavior Z x expected in a doped Mott 
insulator. The lower panel of Fig. [Ml shows that the dop- 
ing dependence of the low frequency optical conductivity 
weight is essentially the same as that of the nodal-sector 
Z. 



VIII. SUMMARY 

In the 8-site DCA approximation to the solution of 
the two dimensional Hubbard model, the doping-driven 
Mott transition occurs in an orbitally selective manner. 
As doping is reduced a pseudogap opens in the region of 
momentum space centered at (0, tt) (and (tt, 0)) while the 
sector containing the zone-diagonal remains ungapped 
down to much lower dopings. In this paper we have 
explored some of the implications of this pseudogap for 
observables. We find that the gap is apparently tied to 
the Fermi level, fills in rather than closes with increas- 
ing temperature, and produces many of the qualitative 
features observed in experiment. 

Fig. [15] plots the doping dependence of the pseudogap 
magnitude as defined from the peak to peak separation in 
the electron spectral function. As is seen in experiment^ 
the pseudogap magnitude is approximately linear in dop- 
ing. The extrapolation to a; = (see Fig. [TSl) gives 
lim2;_>o 2Ap(3(x) = l.lt rather smaller than the a; = 
sector C gap ~ 2.6t characteristic of the half filled insu- 
lator. The pseudogap is thus a new phenomenon, not a 
remnant of the gap occurring at half filling. 

For X > 0.11, a gap is not visible in our calculated 
spectral functions; however, the linear extrapolation of 
the points in Fig. [T5| indicates that for U = 7t studied 
here the critical value at which the gap would close is 
X ^ 0.17. One possible inference is that the temperatures 
T > t/20 accessible to us are too high to enable the gap 
to be seen, and that calculations at lower temperatures 
would reveal a gap in the range 0.11 < a; < 0.17. Indeed, 
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FIG. 14: Upper panel: doping dependence of the quasi- 
particle residue Z calculated for sector B containing the 
(7r/2,7r/2) point and sector C containing the (0, tt) point 
from analytically continued self-energies at inverse tempera- 
ture /3 — 20/t ~ 200 K. The filled symbols represent dopings 
X > 0.08 for sector B and x > 0.18 for sector C where the 
self-energy is Fermi-liquid-like as defined in the text. Open 
symbols are mathematically defined from our data but we be- 
lieve are not physically meaningful because the scattering rate 
is too large. Lower panel: approximate proportionality of op- 
tical spectral weight and Z. Trinagles (black online): integral 
from to a; = 2t of calculated optical conductivities, divided 
by renormalization factor Z of sector B, plotted against dop- 
ing X. Squares (red online): integral from to a; = 2t of 
sector B contribution to conductivity, divided by renormal- 
ization factor Z of sector B and plotted against doping. 



in experiment pseudogap effects are visible at dopings > 
0.11 but only at temperatures < 200 K (see e. g. Ref . fisl. 
Fig. 63 or Ref. Isol) . i.e. less than the lowest temperature 
accessed in this study Extending our results to lower 
temperatures is therefore of interest. It is also interesting 
to note that this larger doping is comparable with the 
value X = 0.15 at which a quantum critical point was 
reported from analysis of the electron lifetime in a model 
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with U = 6t4i 

Fig-HSlstiould however be interpreted with caution. At 
low dopings (see e. g. the x = 0.05 results in Fig. (2) the 
low T limit of the spectral function exhibits a clear gap 
(region of vanishing density of states) and the peak to 
peak distance plotted in Fig. [15] corresponds well to the 
gap edge. As doping is increased the gap fills in and it is 
less clear whether as T — the density of states would 
develop a true gap, or whether for x > 0.1 the peaks 
represent the boundary of a region with reduced, but 
non-vanishing density of states. If the second possibility 
occurs then a; ~ 0.1 is a critical doping separating a low 
doping region where sector C is gapped from an inter- 
mediate region where sector C has a non-zero, but possi- 
bly suppressed Fermi level density of states. This latter 
possibility is suggested by the imaginary axis analysis of 
previous papers^i^ which defined the critical doping for 
the orbital selective transition in terms of the chemical 
potential at which carriers were first added to sector C. 
This chemical potential value implies a critical x ^ 0.1. 

A closely related question concerns the possible forma- 
tion of Fermi arcs. A natural interpretation of the results 
is that as doping is decreased the pscudogap first forms 
at (0, tt) / (tt, 0) and with further decrease of doping an in- 
creasing portion of the Fermi surface is gapped, leaving 
a "Fermi arc" whose width is doping dependent. In this 
interpretation the sector B finding of a scattering rate 
which becomes large as doping is decreased would rep- 
resent an average over a region (increasing as doping is 
decreased) where the Fermi surface is gapped and a Fermi 
arc region (decreasing as doping is decreased) with good 
quasiparticles. Analysis of this possibility requires a finer 
momentum resolution than is presently available to us. 
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FIG. 15: Points; pseudogap size determined from electron 
spectral function as a function of hole doping x for U = 7t 
and inverse temperature /3 = 20/t ~ 200 K. Line: linear fit 
to results. 

While differences remain on the quantitative level be- 
tween calculation and experiment and indeed between 
calculations performed on different clusters, the re- 



sults indicate clearly that the Hubbard model at in- 
termediate correlations and low dopings does exhibit a 
pseudogap with many of the features exhibited by the 
experimentally-defined high energy pseudogap. A tran- 
sition to a phase with long ranged order is not necessary 
to produce the effect. The calculation does not repro- 
duce many of the lower energy anomalies which may be 
associated with onset of significant superconducting, ne- 
matic, or orbital current order, perhaps because these are 
long wavelength effects beyond the scale provided by the 
cluster sizes or because they involve physics beyond the 
one-band model we are presently able to study. 
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Appendix A: Determination of Fermi-liquid regime 
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FIG. 16: Imaginary part of the B-sector self-energy (solid 
line, red on-line) and real part of the _B-sector self-energy 
(dashed line, blue on-line) at dopings x — 0.13 (upper panel) 
and X = 0.08 (lower panel), as a function of frequency. Note 
the different y axis scales. 

In this appendix we provide details of the analysis we 
use to determine whether the system is in a Fermi-liquid 
regime. Representative results for ReE(a;) and Iml](aj) 
are shown in Fig. 1161 In a Fermi-liquid, the real part of 
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the self-energy is linear in frequency (at low frequency), 
and the imaginary part is not too big. To formulate a 
quantitative criterion we first determine the range around 
a; = over which ReE(cj) is linear. This range is bounded 
at the lower end by a frequency lol and at the higher end 
by ujh- For the x = 0.13 data we have wl = —0.08, ljh = 
0.08, for X = 0.08 we have ujl = -0.08, ujh = 0.02. Then 
we compute the change in quasiparticle energy {ujh — 
E{ojh)) - {ujl - (finding 0.56596 for x = 0.13 

and 0.42048 for x = 0.08). Finally we compare the result 
to — 2ImS(a; = 0) which changes considerably between 
dopings. From this comparison we see that sector B at 
X = 0.13 is well within the Fermi-liquid region while x = 
0.08 is on the border. 



Appendix B: Real Part of Self Energies 

For completeness we show in this appendix the self- 
energies not discussed in the text. Fig. [T71 presents the 
imaginary part of the self-energy in sector B contain- 
ing the (±7r/2, ±7r/2) point. The real parts of the self- 
energies in sector B and sector C over an intcrmeidate 
frequency range are shown in Fig. 1181 
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FIG. 18: Real part of the self-energy in sector B (upper panel) 
and C (lower panel) as a function of frequency, for dopings 
indicated, corresponding to the imaginary parts of the self- 
energies displayed in Fig. [13] and Fig. (TT] 



FIG. 17: Imaginary part of the self-energy in sector B as a 
function of frequency, for dopings indicated, corresponding to 
the real part in Fig. 1181 
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